Enhanced Rare Pion Decays from a Model of MeV Dark Matter 
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A model has been proposed in which neutral scalar particles x, of mass 1 — 10 MeV, annihilate 
through the exchange of a light vector boson U, of mass 10 — 100 MeV, to produce the 511 keV line 
observed emanating from the center of the galaxy. The \ interacts weakly with normal matter and 
is a viable dark matter candidate. If the U -boson couples to quarks as well as to electrons, it could 
enhance the branching ratio for the rare decay 7r° — > e + e~ . A recent measurement by the KTeV 
Collaboration lies three standard deviations above a prediction by Dorokhov and Ivanov, and we 
relate this excess to the couplings of the t/-boson. The values are consistent with other constraints 
and considerations. We make some comments on possible improvements in the data. 
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Astronomical observations over the past century have 
shown that approximately 20% of the universe is made 
of dark matter. Since dark matter has only been de- 
tected through its gravitational interactions, several of 
its properties, including the mass of dark matter par- 
ticles and their interactions with Standard Model (SM) 
particles, remain completely unknown. While the WIMP 
(Weakly Interacting Massive Particle) with mass of or- 
der the electroweak scale is arguably the most popular 
candidate at this time, the possibility that dark matter 
could be lighter than 100 MeV has attracted much at- 
tention recently. A model proposed by Boehm et al. [l| 
postulates a neutral scalar dark matter particle x with 
mass 1 — 10 MeV Q, which annihilates to produce elec- 
tron/positron pairs: XX ~ y e + e~ (alternately, x could 
be a fcrmion [3]). The excess positrons produced in 
this annihilation reaction could be responsible for the 
bright 511 keV line emanating from the center of the 
galaxy [j, as more conventional astrophysical explana- 
tions have failed to explain both the intensity and shape 
of this line. Boehm proposes two particles which medi- 
ate x annihilation: a neutral vector boson U, with mass 
mjj ~ 10 — 100 MeV, and a heavy fermion with mass 
> 100 GeV. The U boson is needed to explain the relic 
dark matter density, while the F fermion is necessary to 
account for the observed rate of positron annihilation in 
the galactic center For early incarnations of the U 
boson, see [6|. 

The rare decay 7r° —> e + e" has long posed an interest- 
ing problem in the theory of strong interactions. Since it 
is suppressed at tree-level, the rate is very small, which 
provides an opening for indirect effects of new physics to 
appear. A new precise measurement of B(tt° — > e + e~) by 
the KTeV Collaboration exceeds the most recent theoret- 
ical calculation. In this Letter, we examine the possibility 
that this excess can be explained by the exchange of an 
off-shell U boson of the type proposed in the light dark 
matter model, and we compare the constraints obtained 
on the coupling constants to other constraints. 



The tree-level decay ir° — > 7* — > e + e _ is forbidden be- 
cause the pion is a pseudo-scalar particle and the massless 
photon has no longitudinal component. Consequently, 
the lowest-order SM contribution is a one-loop process 
with a two-photon intermediate state. The suppression 
of the amplitude by a factor of a 2 and by helicity con- 
servation leads to an extremely small decay width, thus 
allowing for even tiny effects of new physics to be de- 
tectable. 

The KTeV-E779 Collaboration recently published a 
new measurement of i?(7r° — > e + e~) They nor- 

malized this branching ratio to the Dalitz decay mode 
B(tt° — > e + e~7), taking both sets of decays from a large 
sample of Kl — > 3ir° events. The similarity of the sig- 
nal and normalization channels serves to minimize accep- 
tance and efficiency uncertainties, and the remaining sys- 
tematic uncertainties are mainly external, deriving from 
the rate of Dalitz decays and the parametrization of the 
pion form factor. The KTeV Collaboration report the 
value 



■) = (7.48 ±0.29 ±0.25) x 10" 



after extrapolating from the selected to the entire kine- 
matic region. The first error is from data statistics alone, 
while the second is the total systematic error. 

The most recent theoretical estimate of the ir° — > e + e~ 
width was completed by Dorokhov and Ivanov, who ob- 
tained 



B SM (7T° 



(6.2 ±0.1) x 10" 



and noted the discrepancy with respect to the KTeV 
measurement @. The dominant theoretical uncertainty 
comes from the hadronic form factors, expressed through 
a dispersion relation in terms of a subtraction constant. 
The authors of [8J] estimate the value of the subtraction 
constant by using an assumed monopole functional form 
and comparing with CLEO data. They further com- 
pare their estimate with results derived from the op- 
erator product expansion, QCD sum rules, generalized 



/„■ is the pion decay constant, and p^ is the ir four- 
momentum, p 2 = m 2 . (See the Appendix n for details). 

To obtain the full amplitude for ir° — ► e + e~, the U 
boson matrix element is combined with the Standard 
Model amplitude for n° — ► e + e _ 8] and summed over 
the outgoing electron and positron spins. The partial 
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is computed from the expression for 



the two-body decay, 
u, d 

7T° J 




(3) 



where \p\ is the t^rge-momentum of one of fehe outgoing 
particles, and is equal to approximately neglecting 
the electron mass. 



III. BOUNDS ON {/-QUARK COUPLINGS 

We interpret the positive difference B mcas — B SM — 
(1.3 ± 0.4) x 1(T 8 as the contribution of Mu in Eq. ©. 
Taking the known pion and electron masses, f n o = 130 ± 
5 MeV and r^o = (84 ± 6) x 10~ 18 s 0, we find 



(9 u A-9 d A )g e , 



1)17 



(4.0 ± 1.8) x 10~ 10 MeV -2 . (4) 
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FIG. 2: Constraints on the couplings as a function of mjj. The 
thick (red) line labeled "7r°" shows our result Eq. Q. The 
dotted line comes from constraints on (g — 2) e and the dashed 
line from v-e scattering |J . The solid (blue) lines labeled with 
"n MeV" are based on relic density calculations 



In order to make contact with other constraints on this 
model, we assume, as an illustration, that the electron 
coupling and the difference in quark couplings are equal, 
i.e., g\ — g\ = g\ = Qa- This choice is arbitrary, but 
one might naturally expect such a relation to hold within 
an order of magnitude; a more precise relation requires 
a specific model for the fermion charges under U(l)u, 
which is beyond the scope of this Letter. With this as- 
sumption, 



g A = 2.0±g;f x 10~ 4 



mu 



10 MeV 



(5) 



where the asymmetric error bars come from taking the 
square root of Eq. (H|) . Fig. [2] shows this constraint as 
a thick line labeled "7r°". If a given model specifies a 
different relation between g\ — g A and g e A , then this line 
will move vertically in the plot. 

Fayet has derived other bounds on the coupling of 
U bosons to quarks and leptons from a variety of pro- 
cesses @, and some of these are shown in Fig. [2] The 
dashed line labeled u (g — 2) e " indicates his constraints on 
the axial coupling of U to electrons derived from mea- 
surements of the anomalous magnetic moment of the 
electron; the region above this line is excluded. Con- 
straints from kaon decays, as well as (g — 2) M [{|, can 
be evaded if we assume that couplings to second and 
third generation fermions are suppressed. Neutrino- 
electron scattering can provide a relatively severe con- 
straint @ , but may be evaded if the coupling to electrons 
is largely right-handed. Finally, the three solid lines la- 
beled "1 MeV," etc., show constraints on the total U — e 
coupling ftot = \J {fy) 2 + (/I) 2 from the dark matter 
relic density [§], assuming C x = 1, for three hypotheti- 
cal values of the x mass. The regions above these lines 
correspond to smaller values of C x . 

The curves in Fig. [2] show that our values for the 
couplings of the [/-boson to light quarks and leptons 
are interesting in the context of the light dark matter 
model, falling in the same order-of-magnitude as other 



constraints. Since M.u depends on a set of coupling 
constants different from the other constraints, the rare 
decay ir° — > e + e~ provides a different view of the phe- 
nomenology of the light U boson. 

We note that the partial width for the electroweak in- 
teraction process 7r° — > Z* — > e + e~ is identical to Eq. © 
with mu replaced by mz and g A d ' e replaced by the anal- 
ogous weak couplings; even with g A d ' e as small 10~ 4 , the 
large mass of the Z boson renders this electroweak pro- 
cess completely negligible compared to the {/-mediated 
decay. 



IV. POSSIBLE IMPROVEMENTS 

The excess depends on one precise measurement and 
one precise calculation, and has a significance of only 3er. 
Improvements to this situation would help sharpen the 
discussion. 

The KTeV measurement depends crucially on the 
branching ratio for Dalitz decays, which was last mea- 
sured in 1981 [ll[ to a precision of 3%. Combined with a 
bubble chamber measurement from 1961 [l2j], the world's 
best value for this basic benchmark decay is precise at the 
level of 2.7% [To| ■ We speculate that modern experiments 
could provide a more precise measurement. For example, 
Na48/2 has reconstructed 9.1 x 10 7 decays K* 1 — > 7r ± 27r° 
as a key part of their pro gram to study CP violations in 
charged kaon systems [13| . Their trigger does not appear 
to have suppressed Dalitz decays. As these decays take 
place in a long vacuum tube, there is no background from 
converted photons. The resolution on photon energies 
and position is excellent, allowing a direct reconstruction 
of the vertex position. Combined with the charged pion, 
the mass resolution is 0.9 MeV/c 2 . Track and photon 
reconstruction efficiencies have been measured directly 
from the data. Perhaps the acceptance ratio for one 
charged track plus five electromagnetic clusters over one 
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charged track plus four electromagnetic clusters could be 
estimated by simulations at the 1% level. If so, a much 
better measurement of B(ir a — » e + e~~f) would appear to 
be possible. Clearly, other experiments might endeavor 
to update this branching ratio measurement. 

Another approach to testing the SM with better pre- 
cision might be to compare the measurement of the ratio 
of branching ratios, B(tt° — ► e + e~)/ B(ir° — > e + e _ 7), to 
the theoretical prediction. Clearly this ratio is experi- 
mentally much more precise, and we expect the theoret- 
ical uncertainties coming from the hadronic form factor 
to be reduced as well. 



Of course, a direct measurement of £>(-7r c 



") by 



another experiment would be extremely interesting. 

We considered possible enhancements to leptonic rj de- 
cays provided by a light [/-boson. Following the same 
calculation as above, we find 



Mr 



3(9 e A(9 U A+9 d A-2g s A )f*) 



[u^^vjpf, 



Assuming that the couplings are of the same order as in- 
dicated in Eq. (JSJ), we find B(rj — * e + e~) ~ 10~ 9 , which is 
of the same order as the improved unitary bound [|[ , and 
much smaller than the experimental bound. 77 mesons 
are heavy enough to decay to muons, and we estimate 
B(r/ —> ~ 2.0 x 10~ 5 . This prediction is nearl y an 

order of magnitude larger than the measured value [14j . 
B(rj -> m + aO = (5-7 ± 0.9) x 10~ 6 . If the [/-boson ex- 
ists, then its couplings to muons must be smaller than its 
couplings to electrons, or the combination of quark axial 
couplings (g% + g d A - 2g s A ) is smaller than (g\ - g d A ). 



APPENDIX: MATRIX ELEMENT 
CALCULATION 



The U -boson exchange couples to both vector and axial 
vector currents of quarks, evaluated between a tt° state 
and the vacuum (see Figure [lj, 

(0| {uY l (<# + 7155a) u + *t (?v + 75<?a) d} \ir°) . _ 
We can reorganize the four terms into combinations with 
definite parity and transformation under strong iso-spin, 
SU(2)v- From there, it is easy to see that the only com- 
bination which contributes is the current proportional to 
the third component of axial iso-spin, 



(0|- {u^ 5 u - d^ lb d) Itt } = Up» . 



where p M = p^ + +P^- is the momentum of the pion. Note 
that the pion state singles out the axial couplings in the 
combination g A — g A . 

The hadronic vertex is contracted through the U bo- 
son, which also carries momentum p into the leptonic 
vertex, 



M = 



(9 u A-9 d A )g e A 



[w7 M 7 5 v] 



9^ 



where u and v are the electron and positron spinors, re- 
spectively. For an on-shell pion, p^p^ = mj, and we 
obtain 



V. SUMMARY 

The 3<T excess of the KTeV measurement of i?(7r° — > 
e + e _ ) over the most recent calculation by Dorokhov and 
Ivanov prompts considerations of a new physics contri- 
bution to this SM-suppressed decay. The light neutral 
vector U boson proposed in the light dark matter model 
provides a good basis for calculating such a contribution. 
We carried out such a calculation and find that couplings 
of the [/-boson to electrons and light quarks should be on 
the order of 2 x 10~ 4 , for my = 10 MeV. Such small cou- 
plings are consistent with the expectations of the light 
dark matter model, and with other constraints coming 
from leptonic processes. 



M= (9 d A-9l)9 A U { - ri5v]p ^ 
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